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ABSTRACT 

Scale  analysis  indicates  that  five  nondimensional  parameters  (Ro1,  t,  y,  X  and  *X)  characterize  the  disturbance 
generated  by  the  steady  flow  of  a  uniform  wind  (Uo,  V<>)  incident  on  a  mountain  ridge  of  width  a  in  an  isothermal, 
uniformly  rotating,  uniformly  stratified,  vertically  semi-infinite  atmosphere.  Here  y  =  ho/ Hr  is  the  ratio  of  the 
mountain  height  ho  to  the  deformation  depth  Hr  =  fa/N  where  /is  the  Coriolis  parameter  and  N  is  the  static 
buoyancy  frequency.  The  parameters  X  =  Hr/H  and  xX  are  the  ratios  of  Hr  to  the  density  scale  height  H  and 
the  potential  temperature  scale  height  H/k  respectively.  There  are  two  Rossby  numbers:  One  based  on  the 
incident  flow  that  is  parallel  to  the  mountain,  e  =  VD /fa,  and  one  normal  to  the  mountain,  Ro  =  UB/fa.  If  Ro2 
4  1,  then  the  mountain-parallel  flow  is  in  approximate  geostrophic  balance  and  the  flow  is  semigeostrophic. 

The  semigeostrophic  case  reduces  to  the  quasi-geostrophic  one  in  the  limit  as  y  and  <  tend  to  zero.  If  the  flow 
is  Boussinesq  (X  =  0),  then  the  semigeostrophic  solutions  expressed  in  a  streamfunction  coordinate  can  be 
derived  from  the  quasi-geostrophic  solutions  in  a  geometric  height  coordinate. 

If  the  flow  is  anelastic  (X  «=  1),  no  direct  correspondence  between  the  two  approximations  was  found.  However 
the  anelastic  effects  are  qualitatively  similar  for  the  two  and  lead  to:  (i)  an  increase  in  the  strength  of  the  mountain 
anticyclone,  (ii)  a  reduction  in  the  extent  (and  possible  elimination)  of  the  zone  of  blocked,  cyclonic  flow,  (iii) 
a  permanent  turning  of  the  flow  proportional  to  the  mass  of  air  displaced  by  the  mountain,  and  (iv)  an  increase 
in  the  ageostrophic  cross-mountain  flow.  The  last  result  implies  an  earlier  breakdown  of  semigeostrophic  theory 
for  anelastic  flow  over  topography. 

Apart  from  a  strengthening  of  the  cold  potential  temperature  anomaly  over  the  mountain,  the  presence  of  a 
finite  potential  temperature  scale  height  (i.e.,  x  nonzero)  does  not  significantly  alter  the  flow  solution. 


1.  Introduction 

The  important  problem  of  airflow  over  mountains 
comprises  a  variety  of  flow  regimes.  Gill  (1982)  pro¬ 
vides  a  clear  review  of  this  literature.  For  mountains 
of  large  horizontal  scale,  the  earth’s  rotation  plays  a 
major  role  in  determining  the  flow  structure.  Meteo¬ 
rological  studies  of  this  regime  have  their  origins  in  the 
work  of  Chamey  and  Eliassen  (1949)  who  presented 
the  first  quasi-geostrophic  analysis  of  airflow  over 
mountains.  Much  of  the  subsequent  work  focused  on 
the  ability  of  the  mountain  anticyclone  to  create  a  stag¬ 
nation  point  and  form  a  Taylor  column  (Bannon,  1980, 
reviews  the  barotropic  literature  on  Taylor  column 
formation.)  Usually1  quasi-geostrophic  theory  treats  the 
lower  boundary  condition  in  linearized  form. 

Robinson  (1960)  introduced  a  formulation  for  ro¬ 
tating  flow  over  topography  that  includes  the  nonlinear 


1  One  exception  is  the  study  of  Buzzi  and  Speranza  (1979). 
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lower  boundary  condition  but  retains  the  important 
quasi-geostrophic  feature  of  the  filtering  of  internal 
gravity  waves.  Jacobs  (1964),  Merkine  (1975),  Merkine 
and  Kalnay  ( 1 976),  Pierrehumbert  (1985)  and  Blumen 
and  Gross  (1986)  have  refined  this  Boussinesq  analysis 
which  represents  an  application  of  semigeostrophic 
theory  (Hoskins,  1975)  to  mountain  flows.  The  present 
study  seeks  to  extend  these  analyses  to  include  anelastic 
effects.  It  represents  a  generalization  of  the  quasi-geo¬ 
strophic  results  of  Smith  ( 1979)  and  Bannon  ( 1986)  to 
finite-amplitude  mountain  ridges. 

A  mountain  immersed  in  an  incident  zonal  wind 
on  a  uniformly  rotating  reference  frame  experiences  a 
lift  force,  L,  due  to  the  pressure  differential  across  the 
mountain  which  tends  to  push  the  mountain  in  a  me¬ 
ridional  direction: 

L  =  -y*  f  pNdA  =  -  f  ^  dV ,  (1.1) 

Js  J  mv  ay 

where  p  is  the  fluid  pressure,  N  is  the  unit  vector  normal 
to  the  elemental  surface  area  dA  directed  toward  the 
fluid,  y  is  the  unit  vector  in  the  meridional  (cross- 
stream)  direction,  S  is  the  mountain  surface  and  MV 
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is  the  mountain  volume.  The  transformation  of  the 
surface  integral  to  an  integral  over  the  volume  of  the 
mountain  utilizes  the  scalar  version  of  the  divergence 
theorem  and  the  Archimedean  concept  of  replacing 
the  solid  body  with  the  equivalent  fluid  (see  Batchelor, 
1967,  pp.  16-17).  If  the  uniform  incident  flow  is  in 
geostrophic  balance  and  has  amplitude  U0,  then  the 
contribution  of  the  ambient  pressure  gradient  to  (1.1) 
is 

L  =  MDX  JV, o,  (1.2) 

where /is  the  constant  Coriolis  parameter  and  MD  is 
the  mass  of  air  displaced  by  the  mountain 

MD  =  I  pdV. 

J  MV 


As  noted  by  Smith  (1979),  a  mountain  induced  pres¬ 
sure  field  could  potentially  generate  an  O(h02)  lift 
(where  ho  is  the  mountain  height)  but  this  cannot  bal¬ 
ance  the  O(ho)  lift  given  by  ( 1 .2).  Invoking  airfoil  the¬ 
ory,  Smith  argued  that  the  lift  ( 1 .2)  generates  a  far-field 
circulation  of  strength 


/X  MD 
MC 


(1.3) 


where  MC  is  the  mass  per  unit  area  of  an  air  column. 
As  (1.3)  indicates,  only  atmospheres  with  finite  MC 
will  possess  a  nonzero  circulation.  For  an  infinitely  long 
mountain  ridge,  a  linear  string  of  circulations  are  gen¬ 
erated;  their  net  effect  is  to  produce  a  meridional  de¬ 
flection  of  the  flow: 


geostrophy.  Section  3  presents  the  closed-form  solution 
of  the  nonlinear  problem  and  the  details  of  its  numer¬ 
ical  solution.  In  addition  the  formula  ( 1 .4)  for  the  per¬ 
manent  turning  is  derived  analytically.  Section  4  de¬ 
scribes  the  results  for  the  standard  anelastic  case  while 
section  5  discusses  the  modified  case.  An  appendix 
compares  the  analytic  quasi-  and  semigeostrophic  so¬ 
lutions  in  the  Boussinesq  case. 

2.  The  model 

The  basic  state  model  atmosphere  consists  of  an  iso¬ 
thermal  compressible  gas  in  hydrostatic  balance.  The 
density/pressure  scale  height  is  H  =  RT/g  where  R  is 
the  ideal  gas  constant,  T  the  uniform  temperature,  and 
g  the  acceleration  due  to  gravity.  The  constant  potential 
temperature  scale  height  is  H/k  =  CpT/g  and  implies 
a  uniform  buoyancy  frequency  N.  Here  Cp  is  the  spe¬ 
cific  heat  capacity  at  constant  pressure  and  k  =  R/Cp. 
A  mountain  ridge  of  width  a  and  height  ho  lies  at  the 
origin  of  a  rotating  Cartesian  coordinate  system  with 
constant  Coriolis  parameter  /  A  uniform  wind  is  in¬ 
cident  on  the  mountain  with  components  Uo  normal 
to  and  V0  parallel  to  the  ridge.  Figure  1  summarizes 
the  model  physics  and  geometry. 


a.  Basic  equations 

The  equations  of  motion  describing  the  inviscid,  hy¬ 
drostatic,  anelastic  flow  are 


Av  =  - 


/X  MDL 
MC 


(1.4) 


d<f> 

dx’ 


(2.1a) 


where  Av  =  v  (downstream)  -v  (upstream)  and  MDL 
denotes  the  displaced  mass  per  unit  length  of  the  ridge. 

Bannon  (1986)  confirmed  the  validity  of  the  far-field 
circulation  for  the  anelastic  quasi-geostrophic  approx¬ 
imation.  In  that  theory  the  lower  boundary  condition 
is  linearized  and  the  displaced  mass  is  approximated 
by 

MD  =  PoMV,  (1.5) 

where  p0  is  the  surface  value  of  the  density.  In  the  pres¬ 
ent  study  the  nonlinear  lower  boundary  condition  is 
incorporated  and  it  is  shown  that  discrepancies  with 
the  linearized  condition  arise  from  ( 1 .5)  being  an  over¬ 
estimate. 

Section  2  presents  the  model  formulation  which 
employs  Robinson’s  (1960)  streamfunction  coordinate 
and  follows  Pierrehumbert  (1985)  in  assuming  a  semi¬ 
infinite  atmosphere.  A  scale  analysis  indicates  that  five 
nondimensional  parameters  describe  the  flow.  It  is 
found  that  the  smallness  of  the  square  of  the  Rossby 
number  based  on  the  amplitude  of  the  flow  normally 
incident  on  the  mountain  provides  justification  of  the 
semigeostrophic  approximation  while  a  mountain 
Richardson  number  measures  the  extent  of  quasi- 


FlG.  1 .  A  schematic  illustration  of  the  model.  A  uniform  wind  is 
incident  on  an  infinite  mountain  ridge  of  height  ho  and  width  a  on 
an  /plane.  The  wind  (U0,  V0)  has  components  normal  (U0)  and  par¬ 
allel  (V0)  to  the  mountain.  The  isothermal  basic  state  atmosphere 
has  density  scale  height  H,  potential  temperature  scale  height  H/k 
and  uniform  buoyancy  frequency  N. 
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Odd 
—  (psu)  +  —  (psv)  +  —  (psw)  =  0,  (2.  Id) 

dx  ay  dz 


D  ,  , 

N2+s 


where 


D  d  d  d  d 

—  =  —  +  «r  +  »-  +  w  —  , 
Dt  dt  dx  dy  dz 

ps  =  p0  exp (-z/H), 


(2. le) 

(2.2a) 

(2.2b) 


N2  =  gd  In djdz  =  Kg/H.  (2.2c) 

The  subscript  s  denotes  the  static  basic  state  fields  and 
59  is  the  dynamic  contribution  to  the  potential  tem¬ 
perature.  In  (2.1)  the  usual  meteorological  convention 
holds  for  the  velocity  field  ( u ,  v,  w )  and  the  geopotential 
0  is  related  to  the  pressure  p  by  0  =  p/pfz).  The  influ¬ 
ence  of  the  mountain  appears  mathematically  through 
the  lower  boundary  condition 


Dh 

w  =  —  at  z  =  h, 
Dt 


(2.3) 


where  h(x)  describes  the  mountain  profile. 

We  write  the  steady  solution  as  the  sum  of  the  in¬ 
cident  flow  and  that  induced  by  the  mountain: 


u  =  U0  +  u'a(x,  z),  (2.4a) 

v  =  V0  +  v'(x,  z),  (2.4b) 

w  =  w’(x,  z),  (2.4c) 

0  =  ~f[U0y  ~  Vox ]  +  4>\x,  z),  (2.4d) 


tf[Uoy  ~  Vo*] 

H 


+  g  T  (x,  z). 


(2.4e) 


The  primed  fields  denote  the  disturbance  generated  by 
the  mountain  ridge  and  by  symmetry  are  assumed  to 
be  independent  of  y.  The  subscript  a  on  the  mountain- 
induced  zonal  wind  component  indicates  that  this  field 
is  an  ageostrophic  one.  We  note  that  the  flow  in  the 
absence  of  the  mountain  is  in  geostrophic  balance  with 
the  uniform  pressure  gradient  force. 

Substitution  of  (2.4)  into  (2.1)  and  (2.3)  yields 


du'q  _  d0' 

dt  JV  dx ’ 


dv' 

dt 


+fu=JU0, 


(2.5a) 

(2.5b) 


d#_ 

dz 


K  ,  5ff 
H*  ~  8  9  ’ 


(2.5c) 


d  (9 

t-  (psu)  +  —  (psw')  =  0,  (2.5d) 

dx  dz 


jt  (n2z  +  gj-)  +  ~  F<>i4)  =  0,  (2.5e) 


where 


d  d  d 

,  dh 

w'  =  ii  —  at  z  =  h(x). 


(2.6) 


V/e  note  that  the  system  (2.5)  and  (2.6)  are  fully  non¬ 
linear. 


b.  Scale  analysis 

A  judicious  scale  analysis  provides  much  insight  into 
the  problem  posed  by  (2.5)  with  (2.6).  At  the  outset  we 
note  that  the  mountain  introduces  zonal  variations  to 
the  flow  through  (2.6).  Thus,  we  choose  a  as  the  char¬ 
acteristic  scale  for  x.  The  choice  of  a  corresponding 
vertical  scale  is  ambiguous  as  five  such  scales  arise  in 
the  problem  (Table  1).  However  the  scales  H  and  H/k 
are  intrinsic  scales  of  the  fluid  but  not  of  the  flow.  Sim¬ 
ilarly  the  mountain  height  may  not  represent  a  typical 
dynamical  scale.  Here  we  choose  HR  rather  than  Ht  as 
the  appropriate  vertical  scale  of  the  flow  as  this  defor¬ 
mation  depth  is  consistent  with  the  characteristic  hor¬ 
izontal  scale  a. 

The  effect  of  the  mountain  is  to  deform  the  isen- 
tropes  horizontally  over  a  distance  a  by  an  amount  A 9 
where 

M~h0%1  =  -N2h0. 
dz  g 

Anticipating  an  approximate  thermal  wind  balance  for 
the  mountain  parallel  flow, 

du'  _  1  d_  (gSff\ 

dz  fdx  \  9S  j  ’ 

the  scale,  VT,  of  the  shear  of  the  meridional  wind  is 

v  =  ^  =  N2fl° 

T  fa9s  fa 

Thus  we  choose  v'  ~  0(  VTHR)  =  Nh0  and,  by  ap¬ 
proximate  geostrophy,  0  ~  O(JNah0). 

The  boundary  condition  (2.6)  introduces  a  scale  for 
the  vertical  field  w'  ~  O(U0h0/a)  provides  u'a  ^  U0.  By 
continuity,  (2.5d),  we  choose  u'a  ~  O(U0h0/HR). 

Based  on  the  above  considerations,  we  introduce  the 
following  scaling 

a:  =  ax ",  z  =  Hrz",  h  =  h0h"(x),  (2.7a) 


Table  1.  Characteristic  vertical  scales. 


h0 

mountain  height 

H  =  RT/g 

density  (pressure)  scale  height 

H/k 

potential  temperature  scale  height 

Hr  —fa/N 

deformation  depth  for  mountain  scale  a 

H,  =  U„/N 

deformation  depth  for  inertial  scale  U0/f 
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u'a  =  UonK,  v  =  Nh0v”,  w'  =  U0  ^  tiw",  (2.7b) 


4>'  =  NfaW,  ^  =  N2h06”, 

Us 

into  (2.5)  and  (2.6)  to  obtain 


where 


(2.7c) 


„  ,  dua  d<t> 

Ro  -~r  ~  v  =  ~  —  , 
dt  dx 

(2.8a) 

dv 

p  —  +  u  =  u, 
dt 

(2.8b) 

dd> 

- —  *X0  =  e, 

dz 

(2.8c) 

d  /A  x  d  x 

—  (pu)  +  p-  (pw) 
ox  dz 

=  0, 

(2.8d) 

—  (z  +  p0)  +  /cXp(u  —  eua) 
dt 

=  0, 

(2.8e) 

d  ,  d  ,  d 

—  =  (U  +  flUa)  ~  h  pW  —  , 

dt  ox  oz 

(2.9a) 

U  =  U  +  pMa, 

(2.9b) 

.  A  .  dh 

w  =  (u  +  pwa)  —  at  z  = 
dx 

ph(z). 

(2.9c) 

In  writing  the  system  (2.8)  with  (2.9)  we  have  dropped 
the  double  primes  and  introduced  u  to  denote  the  scaled 
incident  wind  and  p  for  the  scaled  basic  state  density 
where 

m  =  1,  (2.10a) 

p  =  exp(— Xz).  (2.10b) 

The  five  nondimensional  parameters  (R02,  e,  p,  X  and 
k\)  describing  the  flow  are  summarized  in  Table  2. 


This  differs  from  Merkine’s  (1975)  criterion  that  U02/ 
N2D2  4  1  for  i?0  <  1  where  D  is  the  depth  of  the  fluid. 
Clearly,  as  D  approaches  infinity,  D  is  no  longer  a  char¬ 
acteristic  scale  of  the  flow  and  Merkine’s  criterion  loses 
validity.  Replacement  of  D  with  HR  makes  his  criterion 
equivalent  to  (2.1  la)  and  overcomes  some  of  the  dif¬ 
ficulties  inherent  in  his  nondimensionalization. 

Henceforth  we  assume  that  (2.11)  holds  and  that  the 
flow  is  semigeostrophic.  Then  (2.8a)  is  replaced  by 


In  such  a  case  the  solutions  depend  only  on  e,  p,  X  and 
kX  and  are  not  an  explicit  function  of  the  incident  nor¬ 
mal  windspeed  U0. 

d.  The  Boussinesq  case 

Setting  X  =  0  retrieves  the  Boussinesq  case  and  the 
governing  equations  are  only  a  function  of  p.  This  result 
agrees  with  Pierrehumbert  (1985)  where  p  =  RoFr  in 
his  notation..  He  noted  that  p  is  the  Burger  number 
based  on  the  mountain  slope.  The  present  scaling  in¬ 
dicates  that  p  may  also  be  defined  as 

p  =  Vt/N  =  1/VRi,  (2.12) 

where  Ri  is  the  Richardson  number  based  on  the  ther¬ 
mal  shear  of  the  mountain  parallel  flow.  Equation 
(2.9a)  indicates  that  p  is  a  measure  of  the  advection  by 
the  ageostrophic  flow  while  (2.9c)  suggests  that  p  also 
measures  the  strength  of  the  forcing  by  bottom  topog¬ 
raphy. 

Pierrehumbert  (1985)  noted  a  similarity  between  the 
quasi-  and  semigeostrophic  solutions  for  the  zonal 
wind.  The  Appendix  extends  this  correspondence  to 
the  other  flow  fields. 


c.  The  semigeostrophic  approximation 

Inspection  of  (2.8a)  indicates  that  the  semigeo¬ 
strophic  approximation  of  the  mountain-parallel  wind 
being  in  geostrophic  balance  is  valid  provided 

R02<\,  (2.11a) 

P^O(l).  (2.11b) 

Table  2.  Nondimensional  flow  parameters. 

m  =  h0/HR  =  Nh0/fa 
\  =  Hr/H  =  fa/NH 
k  =  R/Cp 

Ro  =  H,/Hr  =  Uo/fa 
t  =  Vo /fa 


e.  The  quasi-geostrophic  limit 

In  the  limit  as  p  and  e  tend  to  zero,  the  set  (2.8)  and 
(2.9)  with  (2.8a'),  reduces  to 


d<t> 

~V  ~  ~  dx  ’ 

(2.13a) 

dv 

a  +  u"~°- 

(2.13b) 

—  -  k\<p  =  e, 

dz 

(2.13c) 

dua  ,  Id 

~Z~  +  x  —  (Pw)  =  0, 
dx  p  dz 

(2.13d) 

dd 

—  +  kXv  +  w  =  0, 
dt 

(2.13e) 
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where 

£=  -JL 

dt  U  dx  ’ 

w  =  ^j-  at  z  =  0.  (2.14) 

dt 


Solving  (2.13)  for  the  geopotential  yields  the  quasi- 
geostrophic  potential  vorticity  equation 


d  d2(j)  1  3  /  „  30  Y 

dt  dx2  p  dz  \P  dz ) 

with  the  boundary  condition 


=  0, 


(2.15) 


d  /30\ 
dt  \dz) 


dh 

dt 


at  z  =  0. 


(2.16) 


The  Green’s  function  solution  for  the  three-dimen¬ 
sional  problem  appears  in  Bannon  (1986)  while  the 
Appendix  treats  the  two-dimensional  Boussinesq  case. 


)  as  x  -*•  -oo,  v  '  7 
d\ p/dz  =  pu  j  (3.5b) 

The  latter  equation  has  the  solution 


\p  =  -(1/A)e  Xz  as  x  — ►  -oo. 

(3.5c) 

Thus, 

z*  =  — (1/X)  ln(-A^) 

(3.6) 

is  the  height  of  a  streamline  far  upstream  of  the  moun- 

tain. 

With  these  considerations  the  set  (2.8)  and  (2.9)  be- 

comes 

— V  =  —d<t>/dx. 

(3.7a) 

pj(\p,  v)  =  p(u  -  u), 

(3.7b) 

d<t>/dz  —  k\4>  —  8, 

(3.7c) 

z  +  p8  +  K\p(y  +  ev)  =  zw, 

(3.7d) 

pv  =  J(yp,  y). 

(3.7e) 

3.  Analysis 

The  method  of  solution  follows  Robinson  (1960)  by 
introducing  a  streamfunction  coordinate  transforma¬ 
tion  to  form  a  second-order  partial  differential  equation 
for  the  vertical  displacement  of  the  streamline.  In  the 
anelastic  case  treated  here  the  equation  is  an  elliptic 
one  with  nonconstant  coefficients  and  we  obtain  so¬ 
lutions  numerically. 


a.  Derivation  of  the  governing  equation 
Introduction  of  a  mass  streamfunction,  \p. 


p.pw  =  — 


df 

dx’ 


(3.1) 


satisfies  the  continuity  equation  (2.8d)  identically.  A 
convective  derivative  then  becomes,  for  example, 

dv  1 

-  =  (3.2) 

dt  p 

where  J(f,  v)  =  y-(V\p  X  Vt>)  is  the  Jacobian. 

It  is  convenient  to  introduce  the  meridional  dis¬ 
placement  y(x,  z)  of  a  fluid  parcel.  Then  by  definition 


and  the  heat  equation  becomes,  using  (2.8b),  (2.9b) 
and  (3.2), 


—  J[\p,  z  +  p8  +  K\p(y  +  eu)]  =  0.  (3.4) 

P 


with  the  boundary  condition 


p  dx 


didh 
—  T~  at 
dz  dx 


z  =  ph(x). 


(3.8) 


Elimination  of  0  and  6  from  (3.7a,  c  and  d)  yields 
dv 


where 


or 


k\v  =  -  0,  (3.9) 

dz  fx  dx 

0  =  +  €t>),  (3.10a) 

Q  =  z  +  p8,  .  (3.10b) 


is  the  total  nondimensional  potential  temperature  field 
excluding  the  linear  variations  associated  with  the  in¬ 
cident  flow  [see  (2.4e)]. 

We  next  introduce  a  streamfunction  coordinate  sys¬ 
tem  (x,  f)  to  replace  the  Cartesian  system  (x,  z ): 


X  —  x,  \{/ =  \f/(x,  z),  (3.11) 


and  the  streamfunction  height  z  =  z(X,  \p)  is  the  new 
dependent  variable.  Use  of  the  chain  rule  indicates  that 


d  _  d  /3^\  d 

dx  dX  \dxj  dip  ’ 


(3.12a) 


d_  =  !df\  d_ 
dz~\dz)&p’ 

and  the  Jacobian  simplifies  to,  for  example, 


Two  other  useful  relations  are 


(3.12b) 


td\ l/\  dv  .  dv 

(3-12c) 


This  equation  implies  that  z  +  p8  +  K\p(r)  +  ev)  is  a 
function,  zx  say,  of  \p  only.  We  determine  z00(^)  in  the 
usual  manner  by  examining  the  flow  far  upstream  of 
the  obstacle.  There  we  assume  8,  y,  v  and  ua  all  vanish 
and 


—  =  -(d*/dx)/(df/dz), 

d_  =  d_  _  /3z\/#\ 

dx  dX  \dX/\dz/  df 


(3.13a) 

(3.13b) 
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Use  of  the  new  coordinate  system  in  (3.9),  (3.7b)  and 
(3.7e)  yields 

dv  dz  1  (dz  dO  dz  d0\ 

KXVty~»\tydX-dXd$)’  a) 


dv  (* 


(de  dz\  d  r/ae  %  \  dz 
axj  dx  \ax  +  K  pV )  6\p 

<hp[pudi) 


~  k\\V 


-($)■ 


(3.18) 


Infinitely  far  upstream  or  downstream  of  the  ridge,  the 
ageostrophic  motions  vanish.  Then,  from  (3.3),  v  =  drj/ 
dx,  and,  from  (2.8b),  dv/dx  =  0,  and  the  meridional 
wind  is  independent  of  height: 

dv 

—  =  0  at  |x|  =  oo.  (3.19) 


The  variation  of  the  meridional  wind  along  a 
streamline  is,  from  (3.14b), 


dX  g  \u 


(3.20) 


(3.14b) 


V  =  u(drj/dX).  (3.14c) 

The  boundary  condition  (3.8)  with  (3.13a)  becomes 

z  =  fih  at  \]/  =  —  1/X.  (3.15) 

A  second-order  partial  differential  equation  for  z(X, 
i p)  emerges  from  (3.14a  and  b)  upon  cross-differentia¬ 
tion.  The  result  is 


(3.16) 


Integration  along  a  streamline  from  X  =  —  oo  to  +oo 
and  then  across  i p  from  the  surface  (\p  =  —1/X)  to  in¬ 
finity  (\p  =  0)  yields 

An  =  -  f°  dp  r  dx(-  -  l)  (3.21) 
/X  J— 1/X  J— oo  \tl  J 

where  we  have  used  the  fact  that  the  permanent  turning 
An  is  independent  of  i p.  We  use  (3.1)  to  rewrite  the 
integrand  in  the  form 


,^-i) 

#  I 


[Xxp  +  e~Xz] 


(3.22) 


It  can  readily  be  shown  that  this  equation  reduces  to 
the  Boussinesq  problem  if  X  =  0.  The  equation  is  elliptic 
provided 

_pu(dQM)  i 

(d0/dx  +  k\ixv)2  4-  i  ; 

We  note  that  the  numerator  represents  the  static  sta¬ 
bility  of  the  flow  while  the  denominator  is  the  square 
of  a  measure  of  the  thermal  wind  shear  [cf.  (3.9)].  Thus 
the  ellipticity  condition  (3.17)  requires  that  the  effective 
Richardson  number  of  the  flow  in  (X,  \p)  space  be 
greater  than  V*.  Using  (3.10a)  and  (3.14c)  we  rewrite 
(3.17)  as 

(/cX/t)2(^)  (w  -  u)2  4 

and  we  conclude  that  the  standard  anelastic  case  (k 
=  0)  is  always  elliptic.  The  cases  investigated  here  all 
satisfied  (3.17). 

b.  Permanent  turning 

Here  we  derive  the  expression  ( 1 .4)  for  the  net  me¬ 
ridional  displacement  of  the  flow.  First  we  note  that 
this  displacement  is  barotropic.  Use  of  (3.1)  and  (3. 10a) 
into  (3.9)  yields 


using  (2.10b).  Interchanging  the  order  of  integration 
in  (3.21)  and  substituting  (3.22)  enables  the  integration 
over  \p  to  be  performed.  Since  z(\p  =  0)  =  oo  and  z(\p 
=  -1/X)  =  gh,  the  expression  for  An  is 

An  =  —  -  f  dx[  1  -  exp (-Xgh)].  (3.23) 

This  result  is  the  nondimensional  version  of  the  per¬ 
manent  turning  (1.4)  for  an  isothermal  atmosphere. 
Inspection  of  (3.23)  indicates  that  An  vanishes  for  the 
Boussinesq  case  (X  =  0)  and  that  An  =  —X  MV  for  the 
anelastic  quasi-geostrophic  (g  =  0)  case.  Finally  we  note 
that  the  permanent  turning  is  not  a  function  of  k. 

c.  Method  of  solution 

The  coefficients  of  the  elliptic  equation  (3.16)  for  z 
are  nonconstant  and  are,  from  (3.10a),  functions  of  n 
and  r).  We  find  expressions  for  these  variables  using 
(3.14b)  and  (3.14c): 

f  (11-2)^  (3.24.) 

v — CO  fJ,U 


v=  r  (AzA 

v — oo  [Ali 

-£.GK 


(3.24b) 


where 


u  =  1  /(pdz/dt).  (3.24c) 

In  writing  the  definite  integrals  we  have  applied  the 
upstream  condition  that  t;  and  v  vanish  at  X  =  -oo. 
Thus  (3. 16)  and  (3.24)  form  a  coupled  system  of  equa¬ 
tions  for  z. 

We  obtain  numerical  solutions  by  solving  (3.16)  and 
(3.24)  iteratively.  At  each  iteration  we  solve  (3. 16)  using 
standard  techniques  for  fixed  coefficients.  We  use  the 
new  approximate  solution  for  z  to  calculate  revised 
coefficients  from  (3.24).  The  technique  is  convergent 
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after  about  10  iterations.  We  take  a  numerical  solution 
to  the  Boussinesq  case  (X  =  0)  which  has  constant  coef¬ 
ficients  as  the  initial  guess.  1 

We  map  the  infinite  strip  domain  (— oo  <  X  <  +oo_, 
—X-1  <  p  <  0)  into  a  rectangle  (—1  <  X  <  +1,  1  >  p 
>  0)  by  the  transformations 

X  =  tanh(0.2x),  (3.25a) 

P  =  (~\p)l/x.  (3.25b) 

A  uniform  grid  of  101  X  51  points  covers  the  2  X  1 
domain  in  (X,  p)  space.  In  physical  space  the  horizontal 
resolution  is  0.1a  near  the  mountain  and  the  first  in¬ 
terior  points  are  ±  1 1 ,5a  from  the  mountain  peak.  The 
vertical  resolution  near  the  lower  surface  is  0.02 HR  with 
the  highest  interior  point  at  z  =  3.9 HR  from  the  surface. 
We  use  second-order  differencing  and  the  trapezoidal 
rule. 

The  boundary  conditions  for  (3.16)  are 

5  =  0  at  X  =  ±1,  (3.26a) 

5  =  h(X)  at  p=l,  (3.26b) 

8  =  8t  at  p  =  0,  (3.26c) 

where 

5  =  (z  -  ZooVji  (3.27) 

is  the  deviation  of  the  streamfunction  height  from  its 
value  far  upstream.  These  conditions  express  the  phys¬ 
ical  constraints  that  5  vanish  far  upstream  and  down¬ 
stream  of  the  mountain  and  that  the  lower  boundary 
correspond  to  a  streamfunction. 

The  nonhomogeneous  boundary  condition,  (3.26c), 
at  the  top  of  the  domain  requires  some  justification. 
The  constant  8T  is  defined  by 

i  r+t  fnh 

<5r  =  vJ-/  dxl  dz^  ('3'28'* 

and  is  the  constant  displacement  over  the  internal  (-/ 
<  x  <  /)  having  the  same  normalized  mass  per  unit 
length  as  that  displaced  by  the  mountain.  For  our  iso¬ 
thermal  atmosphere,  we  can  perform  the  integration 
over  z  analytically  to  yield 

8 r  =  Yi  J_,  dx^1  ~  exP(~  V*Cx))]/Xm,  (3.29) 

which  has  the  limiting  value 

1  C+I 

8t  =  — ;  dxh(x)  as  (X  or  n)  -*■  0,  (3.30) 

21  J-i 

and  218t  is  the  cross-sectional  area  of  the  mountain. 
Condition  (3.28)  is  therefore  consistent  with  the  an- 
elastic  quasi-geostrophic  (Bannon,  1986)  and  Boussi¬ 
nesq  semigeostrophic  (Pierrehumbert,  1985)  results  for 
which  the  displaced  volume  equals  the  mountain  vol¬ 
ume.  For  a  horizontally  infinite  domain,  l  -*■  oo  and 
8t  -*■  0.  Use  of  this  homogeneous  condition  in  the 


numerical  scheme  yielded  solutions  similar  to  those 
with  a  rigid  lid.  This  suggested  that  despite  the  trans¬ 
formation  (3.25a),  the  effective  horizontal  domain  is 
finite.  Solutions  with  /  =  11.5  (corresponding  to  the 
first  interior  grid  point)  in  (3.30)  showed  excellent 
agreement  with  the  anelastic  solutions  for  the  two  lim¬ 
iting  cases.  All  solutions  presented  here  used  (3.26)  with 
1=  11.5. 


4.  Results  for  the  standard  (k  =  0)  theory 


We  refer  to  the  case  with  k  =  0  as  the  standard  theory 
which  ignores  the  finiteness  of  the  potential  tempera¬ 
ture  scale  height.  In  this  case  the  governing  equation 
(3.16)  with  (3.10a)  and  (3.24c)  reduces  to 


dzx  d2z  d_  7dz^\  1  dz 
dip  dX2  dp  \  dp  )  dp 


(4,1) 


This  result  and  the  boundary  conditions  indicate  that 
the  solution  is  symmetric  about  x  =  0  for  symmetric 
mountains.  In  addition  the  solution  is  independent 
of  e. 

Figures  2-6  display  the  flow  fields  for  the  case  of  a 
Gaussian  mountain 


h  =  exp(-x2),  (4.2) 

with  nondimensional  parameter  settings  X  =  1  and  n 
=  0.4.  Experiments  in  which  these  parameters  are 
changed  indicate  that  the  results  vary  monotonically 
with  both  X  and  tx.  The  choice  displayed  here  is  for 
flow  over  strong  topography  with  significant  anelastic 
effects.  Dimensionally  these  values  roughly  correspond 
to  a  mountain  of  height  h0  3  km  and  width  a  ~  800 
km  in  an  atmosphere  with  scale  height  H  =  8  km,  / 
=  10-4  s-1,  and  N  =  10-2  s-1.  The  semigeostrophic  as¬ 
sumption  will  be  a  good  approximation  (e.g.,  Ro2 
<0.10)  for  typical  atmospheric  wind  speeds  (e.g.  Uo 


Fig.  2.  Horizontal  variation  of  the  vertical  displacement,  5  (in 
units  of  h0),  of  the  streamline  from  its  far-field  value  of  one  defor¬ 
mation  depth  (zj  =  fa/N).  Also  shown  is  the  displacement  of  the 
surface  streamline  which  corresponds  to  the  nondimensional  moun¬ 
tain  height  z  =  h(x).  Here  and  in  Figs.  3  through  6,  fi  =  0.4,  X  =  1, 
and  k  =  0,  and  the  labels  AS,  AQ  and  B  denote  the  anelastic  semi¬ 
geostrophic,  and  anelastic  quasi-geostrophic,  and  the  Boussinesq  cases 
respectively. 
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Fig.  3.  The  cross-mountain  wind  at  the  surface  (in  units  of  U0 )  as  Fig.  5.  The  surface  relative  vorticity  (in  units  of  Nho/a) 

a  function  of  absolute  horizontal  distance.  BS  denotes  the  Boussinesq  as  a  function  of  absolute  horizontal  distance, 

semigeostrophic  case. 


<  25  m  s'1).  For  purposes  of  comparison  the  figures 
also  display  the  anelastic  quasi-geostrophic  (denoted 
AQ)  and  Boussinesq  semigeostrophic  (denoted  BS)  so¬ 
lution.  In  cases  where  the  BS  results  correspond  with 
the  quasi-geostrophic  ones  (see  the  Appendix)  they  are 
denoted  B.  The  new  anelastic  results  are  denoted  AS. 

Figure  2  displays  the  displacement  of  the  zx  =  1 
streamline  which  far  upstream  lies  at  the  elevation  of 
one  Rossby  height  above  the  surface.  Comparison  with 
the  displacement  of  the  surface  streamline  (which  cor¬ 
responds  to  the  mountain  height)  indicates  that  the 
amplitude  of  the  displacement  decreases  with  height 
but  broadens  laterally.  Both  the  B  and  AQ  displace¬ 
ments  conserve  the  mountain  volume  per  unit  length 
given  by  (3.30)  while  the  AS  results  conserves  the  nor¬ 
malized  displaced  mass  (3.29).  While  both  anelastic 
cases  predict  greater  displacements  aloft  than  the 
Boussinesq  case,  the  (linearized)  AQ  result  is  an  over¬ 
prediction. 

The  well-known  enhancement  of  the  cross-mountain 
flow  imposed  by  mass  conservation  is  illustrated  in 
Fig.  3-  Both  the  BS  and  AQ  cases  underestimate  this 
enhancement  and  hence  the  associated  breakdown  of 
semigeostrophic  theory  (Pierrehumbert,  1985).  [We 


Fig.  4.  The  total  mountain-parallel  wind  at  the  surface  (in  units 
of  Nho)  as  a  function  of  distance  downstream  of  the  mountain  top. 
The  fields  are  antisymmetric  about  the  origin.  The  crosses  on  the 
right  ordinate  denote  the  asymptotic  far-field  value  for  each  case. 


note  in  passing  that  the  zonal  windspeed  need  not  be 
infinite  for  the  theory  to  fail.  Large  but  finite  speeds 
would  generate  sufficiently  large  advections  that  their 
neglect  in  (2.8a)  would  no  longer  be  justifiable.]  Farther 
down  the  mountain  slope  the  zonal  wind  is  decelerated 
below  its  upstream  value.  Anelasticity  reduces  this  up¬ 
stream  blocking.  The  zonal  wind  is  symmetric  about 
x  =  0.  Thus,  apart  from  the  regions  of  weak  blocking, 
the  flow  is  divergent  upstream  and  convergent  down¬ 
stream  of  the  mountain  top. 

We  plot  the  mountain-parallel  wind  in  Fig.  4  with 
the  value  of  V0  set  to  produce  an  antisymmetric  field 
with  v(x  =  0)  =  0.  The  permanent  turnings  calculated 
numerically  at  x  =  —11.5  are  At;  =  2  V0/Nh0  =  0.12, 
1 .76  and  1 .53  for  the  cases  B,  AQ  and  AS,  respectively. 
Theoretical  values  for  cases  B  and  AQ  are  At;  =  0  and 
Vir  (=1.77),  respectively.  This  good  agreement  with 
theory  verifies  the  numerical  technique  and  gives  an 
estimate  of  the  numerical  error.  As  noted  in  the  Intro¬ 
duction,  the  overprediction  for  the  AQ  case  with  its 
linearized  boundary  condition  is  given  by  (see  1 .4)  the 
ratio 

PoMV/MD  =  1.15,  (4.3) 

and  agrees  with  the  difference  in  the  anelastic  predic¬ 
tions  for  the  meridional  deflection  Ac.  This  agreement 


Fig.  6.  The  mountain-induced  geopotential  at  the  surface  (in  units 
of  Nhofa)  as  a  function  of  absolute  horizontal  distance. 
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provides  a  posteriori  justification  for  the  boundary 
condition  (3.26). 

The  overprediction  by  quasi-geostrophic  theory  is 
further  reflected  in  Figs.  5  and  6  where  one  sees  that 
the  strength  of  the  AS  mountain  anticyclone  is  reduced 
compared  to  the  AQ  case. 

We  next  note  that  some  of  the  features  of  the  u,  v 
and  f  curves  are  interrelated.  Operating  (3.13b)  on  v 
and  using  (3.14a)  with  X  =  0  and  (3.14b)  yields  an 
expression  for  the  relative  vorticity 


dv  _  u  —  u  1  /<9z\2  u 
dx  nu  h  \dx)  u  ’ 


(4.4) 


ture  anomaly  over  the  mountain  due  to  cold  advection 
of  the  ambient  meridional  gradient  (2.4e). 

Numerical  solutions  for  k  =  2/7,  X  =  1,  and  p  =  0.4 
are  very  similar  to  those  displayed  in  Figs.  2-6  for  k 
=  0  and  are  not  shown.  For  example,  the  maximum 
cross-mountain  windspeed  is  2.65  (2.53),  the  minimum 
zonal  windspeed  is  0.92  (0.93),  the  minimum  relative 
vorticity  is  -1.55  (—1.51),  the  maximum  meridional 
windspeed  is  1.02  (1.01),  and  the  permanent  turning 
is  Av  =  —1.49  (—1.53)  for  k  =  2/7  (0).  These  compari¬ 
sons  suggest  that,  to  within  numerical  error,  the  solu¬ 
tions  are  the  same. 


Since  dz/dx  =  pdh/dx  at  z  =  ph(x),  we  have 


6.  Conclusion 


f  = 


dv 

dx 


z  =  ph(x). 


(4.5) 


We  obtain  insight  into  this  expression  by  taking  the 
quasi-geostrophic  (p.  -*■  0)  limit,  to  obtain 


=  ~ua 


at 


z  =  0. 


(4.6) 


We  note  from  (2.13b)  that  the  quasi-geostrophic  result 
holds  at  all  levels.  Since  the  second  rhs  term  in  (4.5)  is 
O(n),  the  first  term  generally  dominates.  Then  anti- 
cyclonic  (cyclonic)  flow  results  for  u  >  u  («  <  «).  In¬ 
spection  of  Figs.  3  and  5  indicate  that  the  region  of 
blocked  flow  is  typically  cyclonic.  Further  the  over¬ 
shooting  of  the  v  field  seen  in  Fig.  4  requires  cyclonic 
flow.  As  noted  by  Pierrehumbert  (1985),  the  blocking 
depends  on  the  mountain  profile.  Sensitivity  experi¬ 
ments  (not  shown)  for  the  anelastic  case  qualitatively 
confirm  his  findings.  For  example,  solutions  for  a  bell¬ 
shaped  mountain  profile  (which  has  gentler  slopes  than 
the  Gaussian  displayed  here)  exhibit  no  region  of 
blocking,  no  cyclonic  far-field  vorticity  and  the  u-field 
approaches  its  asymptotic  value  monotonically. 


5.  Results  for  the  modified  ( k  =  2/7)  theory 

Here  we  retain  a  finite  potential  temperature  scale 
height  and  set  k  =  2/7.  In  such  a  case  the  governing 
equation  (3.16)  must  be  solved  with  (3.10a),  and  the 
solution  is  a  function  of  e.  In  order  to  obtain  a  sym¬ 
metric  solution  for  z  we  choose  e  =  pAv/2.  As  in  section 
4,  this  choice  makes  the  total  meridional  wind  field 
antisymmetric  about  x  =  0. 

Inspection  of  (3.10a)  indicates  that  the  effect  of  a 
nonzero  k  is  measured  by  the  parameter  grouping 

k\h  =  K^q/H, 


This  study  has  presented  an  investigation  of  anelastic 
flow  over  a  finite-amplitude  mountain  ridge  in  a  ro¬ 
tating,  isothermal  atmosphere.  The  results  indicate  that 
use  of  the  linearized  boundary  condition,  consistent 
with  quasi-geostrophic  theory,  underestimates  the 
strength  of  the  cross-mountain  flow  but  overestimates 
the  mountain  anticyclone  and  the  effect  of  the  lift  force 
in  producing  a  permanent  turning  of  the  flow. 

The  lift  force  (1.1)  may  be  generalized  for  the  beta 
effect  and  nonuniform  incident  flow  to 

L=  f  P(zV(y)U(y,  z)dV,  (6.1) 

J  MV 

provided  the  incident  flow  is  in  geostrophic  balance. 
The  permanent  tuming/far-field  circulation  associated 
with  (6.1)  should  be  modified  accordingly.  However, 
inclusion  of  vertical  shear  in  an  anelastic  flow  and  the 
beta  effect  introduces  a  gradient  to  the  ambient  poten¬ 
tial  vorticity  field.  Thus  Rossby  waves  will  comprise  a 
part  of  the  solution. 

It  is  important  to  note  that  the  lift  (6.1)  depends 
crucially  on  the  presence  of  a  nonzero  incident  wind 
Uo.  Thus  even  though  T  in  (1.3)  and  Av  in  (1.4)  are 
not  explicit  functions  of  Uo,  both  the  far-field  circu¬ 
lation  and  permanent  turning  would  not  be  present  in 
an  atmosphere  at  rest. 

A  final  remark  concerns  the  drag  force  D  acting  to 
push  the  mountain  downstream, 

D  =  -x-  I  pHdA  =  -  I  ^  dV.  (6.2) 
Js  J  mv  dx 

If  the  mountain-parallel  wind  v  is  in  geostrophic  bal¬ 
ance,  then  this  expression  may  be  written  as 

D  =  ~(  p{z)fvdV.  (6.3) 

J  MV 


which  is  the  ratio  of  the  mountain  height  to  the  poten¬ 
tial  temperature  scale  height.  For  the  parameter  settings 
X  =  1  and  p  =  0.4  used  in  the  preceding  section,  k\ h 
=  0.114  and  the  effects  are  small.  In  the  quasi-geo¬ 
strophic  case,  k\h  -►  0,  and  as  noted  by  Bannon  ( 1 986), 
the  modified  results  are  identical  to  the  standard  ones 
except  for  the  strengthening  of  the  potential  tempera- 


In  general  the  drag  (6.3)  does  not  vanish  unless  v  is 
antisymmetric  [assuming  h(x)  is  symmetric].  The  anti¬ 
symmetry  displayed  in  Fig.  4  is  achieved  by  choosing 
the  appropriate  value  of  the  upstream  mountain-par¬ 
allel  flow  V0  (=Nh0  Av/2)  such  that  half  of  the  perma¬ 
nent  turning  occurs  on  the  upstream  side  of  the  moun¬ 
tain  and  half  on  the  downstream  side.  Then  the  drag 
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Table  3.  A  comparison  of  the  quasi-geostrophic  and  semigeostrophic  solutions  for  the  Boussinesq  approximation. 


Quantity 

Semigeostrophic 

Quasi-geostrophic 

Streamline  displacement 

5(X,  p)  =  dkh(k)  exp (ikx  -  \k\p) 

b(x,  z)  =  dkh(k)  exp (ikx  -  \k\z) 

Cross-mountain  wind 

u(X,  P)  =  ii/(l  +  fidS/di p) 

u(x,  z)  =  u(  1  -  pBi/dz ) 

Mountain-parallel  wind 

v(X,  P)  =  f  (3S/dp)dX 

v(x,  z)  =  +f  (d&/dz)dx 

Vertical  motion 

vv(X,  p)  =  u(d&/dX) 

w(x,  z)  =  u(dd/dx) 

Cross-mountain  geopotential  gradient 

Sp  .  . 

—  (X,  z)  =  v 
dx 

Relative  vorticity 

(u\(u\(  db  \2 

ftX,  *)  =  »+  -;  - 

\2j\uj\dXj 

f(x,  z)  =  dv/dx 

vanishes  but  the  mountain  has  an  upstream  influence. 
This  upstream  influence,  however,  is  consistent  with 
the  far-field  circulation  (1.3).  We  further  note  that  for 
this  choice  of  F0,  e  =  gAv/2  and  t  is  not  a  free  param¬ 
eter.  For  the  choice  V0  =  0,  e  =  0,  and  the  u-field  is 
asymmetric  (e.g.,  Blumen  and  Gross,  1986).  As  the 
flow  is  inviscid,  adiabatic,  and  without  wave  motion, 
the  implied  nonzero  drag  must  (Bannon,  1985)  be  the 
result  of  flow  transience.  Here  the  transience  is  asso¬ 
ciated  not  with  the  incident  wind  but  rather  with  the 
starting  vortex  shed  during  the  initial  setup  of  the  flow 
and  forever  adverting  downstream. 
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APPENDIX 

The  Boussinesq  Case 

If  the  flow  is  Boussinesq,  X  =  0  and  the  governing 
equation  (3.16)  reduces  to  Laplace’s  equation  for  the 
streamline  displacement  b.  Pierrehumbert  (1985)  pre¬ 
sented  the  solution  for  a  semi-infinite  atmosphere  in 
terms  of  the  Fourier  transform  h{k)  of  the  mountain 
profile  h{x).  Table  3  summarizes  the  solutions  in  (x, 
\p)  space  for  b  and  related  flow  variables.  The  table  also 
displays  the  corresponding  solutions  in  (x,  z)  space  for 
the  quasi-geostrophic  case. 

Inspection  of  Table  3  indicates  that  the  semigeo¬ 
strophic  solutions  reduce  to  the  quasi-geostrophic  ones 
in  the  limit  as  p.  -*■  0.  Note  also  that  since  \p  =  z  +  fib, 
the  coordinates  (X,  \p)  reduce  to  (x,  z)  in  that  limit. 
Even  for  finite  p  the  two  solutions  possess  a  close  cor¬ 


respondence.  In  particular,  knowledge  of  the  quasi- 
geostrophic  solution  is  sufficient  to  determine  the 
semigeostrophic  one  completely.  For  example,  the  so¬ 
lution  for  b,  v  and  0  (=  —b)  have  the  identical  form. 

Comparison  of  the  two  cases  indicates  that  the 
mountain  anticyclone  has  weaker  pressure  gradients 
in  the  semigeostrophic  cases.  The  vorticity  fields  differ 
in  that  the  anticyclonic  region  is  narrower  and  the  far- 
field  cyclonic  flow  is  stronger  for  the  semigeostrophic 
case.  Despite  the  addition  of  vortex  tilting  and  relative 
vorticity  stretching  to  the  semigeostrophic  vorticity 
dynamics,  the  anticyclonic  vorticity  at  the  mountain 
top  (where  db/dx  vanishes)  is  the  same  for  the  two  cases. 
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